AD  668420 


ARPA  ORDER  NO.  1801 


MEMORANDUM 

RM-6695-ARPA 

MARCH  1068 


OPTICAL  HARMONIC  GENERATION 

AT  A  METAL  SURFACE 


R.  P.  Lutomirski 


D  n 


r 


MAYO  1368 

:  *  |  *  i 


PREPARED  FOR: 

ADVANCED  RESEARCH  PROJECTS  AGENCY 


7U 


nono 


SANTA  MONICA  -  CAllfOftNIA 


f  \  *  •  !,<  -If,--' 

CHAPIN  II  O  \i  %  i 

■  f  *'  !•■>  ,>  .S  i- 

.n  I  -  |  ,  I 


ARPA  ORDER  NO.  188-1 


MEMORANDUM 

RM-5595-ARPA 

MARCH  1968 


OPTICAL  HARMONIC  GENERATION 
AT  A  METAL  SURFACE 

R.  F.  Lutomirski 


This  research  is  supported  by  the  Advanced  Research  Projerls  Agency  tinder  Contract 
No.  DAHC15  67  C  0141.  RAND  Memoranda  are  subject  to  critical  review  procedures 
at  the  research  department  and  corporate  levels.  Views  and  conclusions  expressed  herein 
arc*  nevertheless  the  primary  responsibility  of  the  author,  and  should  not  Ite  interpreted 
as  representing  the  oflieial  opinion  or  policy  of  ARPA  or  of  The  RAND  Corporation. 

DISTRIBUTION  STATEMENT 

Distribution  of  this  document  is  unlimited. 


nu 


Mill) 


tit* 


i  >00  is*  in  M 


ft*. 


PREFACE 


This  Memorandum  was  prepared  for  the  Advanced  Research  Projects 
Agency's  high  energy  laser  program.  It  examines  the  interaction  of 
a  high  intensity  laser  and  a  metal ,  and  discusses  the  heating  and 
second  harmonic  generation.  The  discussion  should  be  of  interest 
to  those  concerned  with  high  energy  laser  devices  and  with  the 
propagation  of  intense  laser  beams  through  bounded  media. 
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SUMMARY 

This  Memorandum  presents  a  calculation  of  the  classical  interaction 

t 

between  laser  light  and  a  metal.  Using  a  simple  model  for  the 
statistics  of  the  conduction  electrons  and  their  interaction  with  the 
surface.  Maxwell's  equations  and  the  Boltanann  equation  are  solved 
self-consistently  for  the  fields  in  the  metal,  yielding  the  usual 
Fresnel  solutions  plus  correction  terms  which,  at  optical  frequencies 
and  above  liquid  helium  temperatures,  are  shown  to  be  of  the  order  of  one 
percent. 

When  laser  light  of  sufficient  intensity  is  used,  a  nonlinear 
polarization  is  induced  and  waves  of  twice  the  fundamental  frequency 
are  produced.  For  the  indicated  model  and  for  normal  incidence,  the 
second  harmonic  waves  in  the  metal  are  calculated  and  the  magnitude 
of  the  relative  heating  is  computed.  The  method  for  solving  the 
general  problem  for  an  arbitrary  incidence  angle  is  outlined  and  the 
ratio  of  the  average  energy  flux  reflected  in  the  second  harmonic  to 
the  incident  flux  is  estimated.  The  relevance  to  a  high-intensity 
laser  experiment  is  considered  and  an  error  in  the  computation  of  the 
reflection  coefficient  as  obtained  by  Jha  *  J  is  discussed. 
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I.  INTRODUCTION 


When  a  la.'-«;r  beam  of  sufficient  strength  interacts  with  a  metal, 

a  nonlinear  polarization  is  induced  in  the  metal  and  waves  of  twice 

the  frequency  of  the  incident  radiation  are  produced.  The  amplitude 

of  these  waves  decreases  rapidly  with  distance  in  the  direction  of 

propagation  and,  except  for  very  thin  foils,  there  is  no  transmission 

of  waves  through  the  metal.  However,  for  oblique  incidence,  second 

harmonic  waves  will  be  reflected  from  the  surface.  While  it  is 

difficult  to  perform  high- intensity  experiments  without  damaging  the 

metal,  claims  of  detecting  a  frequency  doubling  of  the  reflected 

(3) 

light  have  been  made. 

A  number  of  calculations  of  the  nonlinear  conductivity  tensor 
have  been  made  based  on  a  free-electron  gas  model.  7  In  the  usual 
theory  it  is  assumed  that  the  current  density  at  any  point  in  the  metal 
la  determined  entirely  by  the  value  of  the  electric  field  at  that  point. 
With  this  assumption  Maxwell's  equations  have  solutions  representing 
simple  harmonic  plane  waves  which  are  exponentially  damped  in  the 
direction  of  propagation.^  These  expressions  for  both  the  fundamental 
and  second  harmonic  solutions  will  be  referred  to  as  the  classical 
results. 

The  basic  assumption  of  the  above  theory  is  not  valid  when  the 
electric  field  varies  appreciably  over  distances  of  the  order  of  the 
conduction  electron  free  path,  i.e. ,  when  the  free  path  is  of  the  order 
of  the  penetration  depth  of  the  electric  field.  In  general  the 
expression  for  the  current  density  takes  the  form  of  a  definite 
integral  involving  the  electric  field  at  all  points  in  the  metal,  and 


Maxwell's  equations  lead  to  an  integro-differential  equation  from 
which  the  electric  field  has  to  be  determined. 

For  the  fundamental  wave  the  present  theory  predicts  a  propagating 

disturbance  in  the  metal  which  is  characterized  essentially  by  the 

classical  index  of  refraction  with  a  small  correction,  plus  a 

"transient"  term  (from  a  branch  cut)  which  is  of  small  amplitude  and 

decreases  rapidly  with  distance  from  the  surface.  However,  for  the 

second  harmonic  the  theory  predicts,  in  addition  to  the  propagating 

second  harmonic  component  of  classical  theory  (with  a  small  correction) 

and  a  small  transient,  a  propagating  disturbance  approximately  equal 

in  magnitude  and  opposite  in  phase  to  the  classical  result  at  the 

surface,  but  more  rapidly  attenuated.  Hence  the  new  term  is  important 

in  computing  surface  effeccs.  The  magnitude  of  these  terras  at  the 

surface  is  shown  to  be  approximately  eE^/n.caj^  «  (4.2  x  10  )  (inks 

units),  where  e  and  m  are  the  electronic  charge  and  mass,  m  is  the 

P 

2  %  16 

plasma  frequency  (ne  /m«o)^  **  1.4  X  10  rad/sec,  c  is  the  vacuum 
speed  of  light,  and  E^  is  the  amplitude  of  the  incident  electric 
field. 

For  the  present  treatment,  the  metal  is  assumed  to  be  semi- 
infinite  in  extent  and  the  problem  is  formulated  as  a  boundary  value 
problem.  The  relevant  equations  and  approximations  are  discussed  in 
Section  II.  In  the  third  section  solutions  for  the  fundamental  wave 
(frequency  tu)  are  found  using  Fourier  transforms,  and  the  inversion 
involves  a  contribution  from  a  pole  and  a  branch  cut.  The  contribu¬ 
tion  from  the  pole  approaches  the  classical  result  as  v^/c  -•  0 
(where  vf  is  the  Fermi  velocity),  with  the  correction  term  a  function 
of  both  frequency  and  the  relaxation  time  of  the  conduction  electrons. 
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This  correction  to  the  classical  refractive  index  is  small  except 

at  liquid  helium  temperatures;  for  light  in  the  infrared  shining  on 

silver  at  room  temperature  Che  correction  is  of  the  order  of  one 

percent.  The  contribution  fro,*  the  branch  cut  is  also  discussed  and 

its  magnitude  is  shown  to  be  a  small  fraction  of  the  pole  contribution 

3  -1 

at  the  surface  (e.g,,  less  than  one  percent  for  light  at  10  cm 
shining  on  silver)  and  to  attenuate  rapidly  with  distance.  Finally, 
the  correction  to  the  linear  reflection  coefficient  is  considered, 
and  is  shown  to  be  negligibly  small. 

The  second  harmonic  problem  for  normal  incidence  is  considered 
in  Section  IV,  and  an  analysis  similar  to  that  of  Section  III  yields 
the  contributions  from  the  poles  and  the  branch  cut.  It  is  demonstrated 
that  for  normal  incidence  no  second  harmonic  components  are  reflected 
from  the  surface.  In  Section  V  the  second  harmonic  problem  for  oblique 
incidence  is  discussed  and  the  nonlinear  reflection  coefficient  is 
estimated.  The  relevance  to  a  high-intensity  laser  experiment  is 
also  discussed,  and,  in  the  light  of  the  present  analysis,  results 
obtained  by  Jha^  are  demonstrated  to  be  in  error. 


II.  PHYSICAL  MODEL 


For  the  purpose  of  this  analysis  the  penetration  depth  of  the 
field  is  assumed  to  be  small  compared  with  the  linear  dimensions  of 
the  specimen,  so  that  it  is  permissible  to  regard  the  metal  as  occupy¬ 
ing  the  half-space  z  i  0.  The  physical  properties  of  the  metal  can 
be  adequately  described  in  terms  of  the  free  electron  model,  according 
to  which  the  valence  electrons  are  able  to  move  about  freely  through 
the  volume  of  the  specimen.  In  the  absence  of  the  light  wave  the 
zero-order  energy  distribution  of  the  electrons  is  taken  to  be  that 
of  a  Fermi  gas  at  absolute  zero.^  Under  the  combined  action  of  the 
applied  electromagnetic  field  and  the  collisions  of  the  electrons  with 
the  lattice,  a  steady  state  is  set  up,  and  the  distribution  function  is 
determined  by  the  Boltzmann  equation.  All  of  the  details  of  the 
collision  processes  are  summarized  for  the  present  purpose  by  specifying 
the  relaxation  time  T,  or,  equivalently,  the  collision  frequency  v  =  T 

The  collisions  are  caused  by  thermal  or  structural  imperfections  in  the 

14  -1 

lattice;  for  example,  in  silver  v  is  of  the  order  10  sec  .  The 

amplitude  of  the  incident  light  (assumed  monochromatic)  is  then  treated 

as  a  perturbation  parameter  in  ''he  Boltzmann  equation. 

It  is  further  necessary  to  make  an  assumption  concerning  the 

reflection  of  the  electrons  from  the  surface.  The  surface  potential 

varies  from  a  value  of  the  order  of  the  Fermi  energy  at  the  surface 

-8 

to  zero  at  a  distance  d  w  10  cm  in  from  the  surface.  The  average 
time  taken  by  an  electron  to  enter  and  leave  the  barrier  region  is 
then  of  the  order  d/v^.  Hence,  provided  tu  «  v^/d  =-  10^  sec  * , 
it  is  permissible  to  replace  the  surface  potential  barrier  by  a  step 


potential--!. e. ,  to  assume  the  conduction  electrons  are  specularly 
reflected  from  the  surface. 

In  discussing  the  anomalous  skin  effect  in  metals,  Reuter  and 

/ON 

Sondheimer'  '  have  shown  their  results  to  be  relatively  insensitive 
to  the  fraction  of  electrons  assumed  specularly  reflected  from  the 
surface,  and  reasonable  results  can  therefore  be  expected  at  higher 
frequencies.  The  mathematical  formulation  is  presented  in  the  next 


section 
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III.  THE  LINEAR  PROBLEM 

The  electron  distribution  f(r,v,t)  satisfies  the  Boltanann 
(9) 

equation  ' 

if  eE  if 

f +(v.v)f  -r  •§--»<«- V  (i) 

Let 

f  -  f  +  S  f  e“iqa3t  (2) 

0  q  q  W 

where  is  the  distribution  of  a  degenerate  Fermi  gas  at  0°K.  The 
fundamental  electric  field  in  the  metal  may  be  written  as 

E  -  E^e'1®*  (3) 

and  is  proportional  to  the  amplitude  of  the  Incident  electric  field  E^. 
Treating  E^  as  a  small  parameter  in  a  perturbation  expansion, 

£1  *  °<Et>'  £2  -  °<E1> 

If  we  linearize  Eq.  (1)  and  seek  solutions  which  vary  sinusoidally  with 
time  and  depend  only  on  the  coordinate  z,  we  obtain 


ifl  .  kii  f  „  eELzlifo 

3z  v  1  rav  ?>v 


where  asr  ■  oo  +  iu. 

The  formal  solution  of  Eq.  (4)  may  be  written  as 


(4) 


df  1<d'z/v  _  z  .  / 

*!<».*>  ■  ST  aT  e  ‘  [f0  E<«> e  1  d«  +  A<s>]  <5> 


Z  X 


loj#z/v 


Then  e 


increases  indefinitely  with  z  for  v  <0.  In  order  that 

2 


f(«°,y)  remain  finite,  f^  must  be  written  as 


fj(z»v)  -  - 


df  io'z/v  00 

^  •  zlz  ■«>< 


-iffl'l/v. 


mv  9v 

2  X 


for  v  <0 
z 


For  v  >0,  write 
2 


«  of  Ua'z/v  _  z 

f,  ■  7—^  e  z  Tf  E(6ie 

1  mvz  dvx  LJ0 


-  -ifl»'|/v 

[j  E(|)e  zd|  +  A(v)j 


The  specular  reflection  condition  at  z  »  0  implies  that 


ft(0,v  »v  »v.)  "  fi  (°»vv  »vv> 

X  x  y  z  *  /  * 


Substituting  Eqs.  (6)  and  (7)  into  Eq.  (8)  yields 


“  -iw'i/v 

A(v)  «  j  E({>  e  d$ 

J0 


and  therefore 


^  df  io'z/v  z  -iaj't/v 

£-^^e  [JoE<0  e 


Icd'I/v 


d5  +  JoE(|)  e  *d|] 


The  current  density  is  given  by 


(6) 


(7) 


(8) 


(9) 


& 
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j  -  -e  J  V 


f  d  v 


(10) 


If  the  functions  P^(ii)  are  defined  by 


*  ®  *®  rv.  3£  iffl'u/v 

Pi(u>  ■  J  dvx  J  dvy  Nr  sT  e 

v  -«  M  -QO  J  *  Q  X  J 


(11) 


then  Pj:(P')  vanishes  unless  i  ■  x  because  f  *  0  at  v^  «  ±  <».  Dropping 
the  subscript  x,  note  that 


o»  *  oo  v  3f  io/u/v 

j  dTx  J  dvy  J  dv2[r  if  e  T  ■  P(ll>  •  P(‘u)  <12> 

*  -oo  *  -oo  J  v  -oo  Z  X  J 


The  total  current  density  can  then  be  written 


2  j  oo  oo 

jx-  -  “  [  JQ  E(0  P(**5)d|  +  J  E(!)  P(^l)d|  +  J  E(|)  P(e-,)d|1 


*  -  ® 


-  ~  [  J  E(5)  P(|z-5|)d|  +  E(|)  P(*f|)df| 


(13) 


The  electromagnetic  fields  are  related  to  the  current  through  Maxwell's 
equations 


and 


V  X  E  -  io£ 


icaE 

V  X  B  -  V.  J - r 

°"  c2 


(14) 


(15) 


Because  of  the  convolution  form  of  Eq.  (13),  it  is  desirable  to  solve  the 
system  by  means  of  Fourier  transforms.  Equation  (13)  is  valid  only  in 
the  half-space  z  £  0,  however,  and  while  several  techniques  exist  for 
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treating  this  type  of  problem,  perhaps  the  simplest  is  a  mathematical 
artifice  first  used  by  Shafranov. ^  The  technique,  as  applied  to 
this  problem,  is  to  assume  that  the  metal  occupies  all  space,  i.e., 
that  Eqs.  (13),  (14)  and  (15)  hold  for  z  <  0  as  well  as  for  z  >  0, 
and  to  introduce  a  Dirac  delta  function  current  source  at  z  ■  0  to 
yield 

VxB+^E»  iiQ[j  +  J6(z)ex]  (16) 

c  u 

Equations  (3),  (14)  and  (16)  yield 

*2  +  ~~2~  +  *  -ijnu  J6(z)  (17) 

The  general  solution  of  Eq.  (17)  is  a  valid  solution  in  the 
half-space  z  >  0  of  the  original  problem.  The  constant  J  must  then  be 
chosen  such  that  the  fields  connect  properly  to  the  vacuum  field  at 
z  =  0  ,  so  that  the  tangential  electric  and  magnetic  fields  are 
continuous  at  the  boundary.  The  result  of  this  procedure  is  a  solu¬ 
tion  which  satisfies  both  the  correct  differential  equation  for 
z  >  0  and  the  proper  boundary  conditions  at  the  surface  of  the  metal. 

Then,  substituting  Eq.  (13)  into  Eq.  (17)  yields 


2  2  • 

“T  +  ^E  -  f  K(z-5)  E(Od?  -  f  K(z+$)  E(£)d{  -  -i^  J6(z)  (18) 

bz2  c  JO  Jo  ° 


where 
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ifl)  m  ®  *  ®  v  <Jf  |u]/v 

K(u) - §-  J  dv  J  dv  J  dvs  -  ^  * 

o 

2 

and  0)p  *  (n^e  /®«Q)  ,  the  plasma  frequency  of  the  electron  gas.  From 
the  symmetry  of  the  "new"  problem,  E(z)  =  E(-z);  hence  Eq.  (18)  may 
be  written  as 


(19) 


-2  2® 

“-§  +  “j  E  -  f  K(z-|)  E(|)d5  -  -iowoJ6(z) 

„  dz  c 

Equation  (20)  can  now  be  solved  by  a  Fourier  transformation. 
Defining  the  Fourier  transform  of  the  function  f(z)  by 


?(k)  *  ^  f(z)e  dz 


(20) 


(21) 


the  transformed  equation  becomes 


[-k2+  («u/c) 2  -  K(k)  ]  E(k)  «  -iauoJ 


(22) 


where 


K(k) 


m 

-  f  K(Z) 


ikz  .  jg  % 
e  dz  -  » 

n  c 
o 


Liy*  Ldvy  Ldv*  [' 


df 


dv  v  -o>7k 

X  z 


]  (23a) 


2 

CD 

flL.  P-.  ■ 

'  2 
n  c 
o 


CD 


(23b) 


and  where  Eq.  (23b)  follows  from  integrating  Eq.  (23a)  by  parts. 

By  the  Fourier  inversion  theorem,  the  inverse  transform  of 
the  function  1!(k)  defined  by  Eq.  (21)  is 


11 


‘o-sj  lw 

W  -Oft 


ikz  . 
e  dk 


end  therefore  the  solution  of  Eq.  (20)  is 


E(z)  - 


lCDiJL.  J  ® 

_o  p 


ikz .. 
e  dk 


2TT  J—  k2-(m/c)2+K(k) 


where  the  integration,  as  in  Eq.  (24),  is  to  be  carried  out  along  the 

real  axis  in  the  k  plane.  The  usual  procedure  for  evaluating  such  an 

integral  is  to  close  the  contour  in  the  upper  half-plane  and  to  apply 

the  theory  of  residues.  Now  assume  f  will  vanish  for  v  greater  than 

o  z 

some  value  v^.  Equation  (23)  then  reveals  that  K(k)  is  not  defined 
whenever  kv^/co'  is  both  real  and  has  magnitude  greater  than  unity, 
i.e. ,  whenever  arg  k  ■  ±  6  =  ±  arg  eo'  ■  ±  tan  *(v/a>)  and  Jk|  ^  Q/v^ 
where  Cl  •  |o/|  ■  (eo +u  )  .  These  conditions  define  the  branch  cuts 
around  which  the  contour  must  l’  deformed.  The  deformed  contour  and 
the  branch  lines  are  indicated  in  Fig.  1,  where  E(z)  is  recovered 
in  the  limit  as  R  -•  co. 

The  integral  over  the  infinite  semicircle  is  easily  shown  to 
vanish.  Then  the  residue  theorem  applied  to  the  contour  yields 


imi  J 
o 


ikz,, 
e  dk 


o  r  e  dk 
2TT  *Br  k2-(m/c)2+K(k) 


[IKZ  _ 

1 — —y-  ] 

j  k -(o>/cr+K(k)J. 


where  each  k^  is  an  Isolated  pole  of  the  Integrand  of  Eq.  (25)  in  the 
upper  half-plane,  the  integral  is  along  the  branch  cut  and  denoted 
by  Br,  and  Res  is  the  residue. 
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CONTRIBUTION  TO  THE  ELECTRIC  FIELD  FROM  THE  POLE 

In  this  section  the  contribution  from  the  residues  to  the  electric 
field  are  calculated.  To  locate  the  poles  in  the  k  plane,  the  roots 
of  the  dispersion  relation 


D(k)  =  k2-(to/c)2+K(k)  *  0 


(27) 


must  be  found.  These  roots  are  discussed  in  Appendix  A.  In  the  limit 

x  »  v  it  is  shown  that  there  always  exists  one  and  only  one  zero  of 

this  expression,  provided  f  is  a  symmetric  function  of  v  ,  and  it 
,  o  z 

kv 

occurs  when  | — ^-7 1  «  1.  The  integrand  of  Eq.  (23)  is  an  analytic 
function  in  the  "cut"  plane  and  lt(k)  can  be  evaluated  analytically 
for  in  »  v  by  expanding  the  denominator  to  yield: 


K(k) 


2 

CD 

01_  _E_ 

car  2 
c 


.2  2 

k  v 


(28) 


The  index  of  refraction  n  =  (kc/u))  is  then  given  by  the  positive 
square  root  of 


2 

n  = 


u,2 

1  "  aim7 


2  2 
CD  CD  V 

1  +-2  _5 
a,'3  c2 


(29) 


Equation  (29)  describes  the  wave  propagation  characteristics  of  a 

metal  as  a  function  of  frequency,  including  the  effect  of  a  finite 

velocity  spread  in  the  distribution  function  for  the  free  electrons. 
~2  2 

As  the  ratio  v^/c  tends  to  zero,  the  above  expression  continuously 
approaches  the  usual  formula  for  a  "cold"  electron  gas 
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n? !  -  ft)  -  1  - 
c  \  <a  J 


cd(otHu) 


(30) 


The  magnitude  of  the  thermal  correction  in  n  ,  4 


2  2 
n  -n 
— — £ 

1  >  18 


(31) 


The  maximum  value  of  the  correction  occurs  at  ai  »  u/  •/!.  Substituting 
in  Eq.  (30)  reveals  that 


L  * 


,3/2 


2  2 
,00  ,  v 

(?)  i 


(32) 


Hence  for  cd  »  u,  only  the  first  few  moments  of  f  need  be  estimated 

to  compute  the  dispersion  relation.  For  lower  frequencies  it  is 

necessary  to  approximate  the  shape  of  the  Fermi  surface  by  an  analytic 

expression  or  to  numerically  integrate  the  expression  for  ST(l:) .  In 

general,  the  energy  surface  can  be  complicated,  and  its  shape  for 

certain  metals  has  been  estimated  from  measurements  of  various  optical 

constants,  soft  x-ray  emission,  Knight  shift,  and  various  transport 

properties. For  rough  order-of-magnitude  calculations  we  can 
~2  2 

aastme  is  0(vf)  where  vf  is  calculated  from  the  free  electron  model 

and  depends  only  upon  the  valence  electron  density.  A  spherical 

energy  surface  of  the  free  electron  model  is  expected  to  be  a  good 

approximation  in  a  number  of  cases,  including  sodium,  silver,  and 

(12) 

certain  alloys  of  copper.  '  While  it  is  not  necessary  to  assume  a 
spherical  distribution  to  proceed  with  parts  of  the  analysis,  if 
some  of  the  integrations  are  performed  the  formulas  become  simpler  and 
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permit  us  to  readily  estimate  the  magnitudes  of  the  various  terms. 
Then  with 


where 


f 

o 


n 

o 


(4/3)tt  v3f 


H(vf-v) 


H(x)  -  1,  x  s  0 


(33) 


-  0,  x  <  0 

2  12  o  8 

we  have  v*  •  j  v*.  For  silver  at  0  C,  vf  «  1.4  x  10  cm/sec, 

a>p  as  1.4  X  10*6  rad/sec,  and  v  «  .5  x  10^  sec  ^  Then  with 

c  «  3  y  10^  cm/sec,  the  maximum  value  of  A  is  .14.  For  poorer 

conductors,  A  will  be  correspondingly  less, 
max 

The  root  of  Eq.  (27)  can  then  be  written  approximately  as 


k.-kc 


(34) 


2  2 

where  cp  »  tan  3tU  fc  is  the  classical  propagation  vector, 

^(to  -3u  )J  c  3-1 

and  d/2  £  .07.  For  example,  for  light  at  10  cm  shining  on  silver, 
q  -  2  X  lO^sec  *,  and  from  Eq.  (31),  d/2  »  .01. 

Returning  to  Eq.  (26) ,  the  contribution  to  the  electric  field 
from  the  pole,  E^,  becomes 


E  «  -  on  J  x  Res 
P  o 


.  lks_ 

[too] 


-o*i  J  e 
_ 2 _ 


ik  z 
o 


<D  J 


(35) 


The  magnetic  field  associated  with  E^  is  given  by 
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1 

B  ( z)  *  7—  - —  e 

-p  Ud  ftz  -y 


Me 

o 


lk  z 

o 


'[■'iff1)  (3)1 


(36) 


The  constant  J  cannot  be  evaluated  until  the  contribution  from  the 
branch  cut  is  calculated.  However,  in  most  practical  cases,  the 
propagation  characteristics  of  (and  Bp)  will  be  shown  to  be  given 
approximately  by  the  classical  electron  theory  (i.e.,  kQ  «  kc) ,  with  a 
correction  of  at  most  a  few  percent  above  liquid  helium  temperatures. 

CONTRIBUTION  TO  THE  ELECTRIC  FIELD  FROM  THE  BRANCH  CUT 


In  the  integration  along  the  branch  cut  in  Eq.  (26) ,  the 
contribution  from  the  small  semicircle  is  easily  seen  to  approach 
zero  as  e  -»  0.  Then  the  contribution  from  the  branch  cut  is 


ikz 

e  dk 


J 

^  2lT  JCl+c2  k2-(oi/  c)2+K(k) 


(37) 


Along  c^,  i.e.,  on  top  of  the  branch  cut,  k  p  e*^  e*e.  Then 
ocVk  «  (OeiS)/(p  e^e**)  -  -  e  ie.  Therefore,  the  pole  in  the  integrand 
in  Eq.  (23)  lies  just  below  the  real  axis  in  the  complex  v^  plane. 
However,  the  pole  can  be  considered  to  lie  on  the  real  ax.  and  the 
path  of  integration  can  then  be  deformed  to  lie  just  above  the  pole. 

The  new  path  will  be  denoted  by  c+,  and  the  integral  so  defined  by 
l?+(p).  Similarly,  along  c  the  path  in  the  plane  can  be  deformed 
to  c  (see  Fig.  2) ,  and  the  corresponding  integral  will  be  denoted  by 
IT  (p) .  Then 


Fig. 2— -Paths  of  integration  for  k+(k)  and  k  (k) 
in  the  second  harmonic 
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2ttE_  ®  izpe*^  i5  .  vf  izpe'-''  i6  . 

_ Br  „  f  _ e  H  e  dp  r  e  e  dp 

iaii  J  Jn/  2  2i6  2,2  “+,  N  J  2  2i6  2,  2  ,  ~- 

o  *C)/v,  p  e  -  to  / c  +  K  (p)  "«  p  e  -  03  /c  +  K 


0/v,  .  16 

f  izpe  i6 

e  dp 


(P> 


. -  ®  .  i6 

io  r  izpe 


•  f,  e 


Jn/vf 


r  2  216  2,2,  «s+,  .nr  2  216  2,2^  .n 

[p  e  -  oi  / c  +  K  (p)jip  e  -  co  /c  +  K  (p)J 


(38) 


where 


03  ®  ®  _9£  v 

^(p)  -  ®  -*2  f  d\J  %  J  d\[jf  ^"n] 

n  c  •'  -®  “  -®  "c,  L  xv - 


(39) 


z  p 


(ir*(p)  -  K~l'(p)  3  can  als°  be  expressed  by  Eq.  (39),  except  that  the 
integration  is  now  performed  along  a  closed  contour  enclosing  the 
pole  (in  the  positive  direction).  Consequently,  according  to  the 
theory  of  residues, 

2 


K  -  K+-  2tri“ 
P 


03  *  „®  ,  df 

^2  J  dvx  f  dvy  (Vx  *r) 
n0c  x  v  -n/p 


(40) 


From  Eq.  (33), 


df 
_ c 

dv 


x 


2  2.2 
v  +v  +v 
x  y  z 


3fwv2+v2-fv2  -  V  ) 

'  v  x  y  z  f/ 


(41) 


2  2 

Then,  changing  to  cylindrical  coordinates  and  substituting  v  -* 

+  (n/p)2  yields 


4- 

K  -  K  - 


- 1  "i  ®  Cv2c-(n/S)2]  B[v2-<n/p)25 


2  p  2  3  f 
c  vf 


(42) 
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Therefore 


JBr 


2  2 

1 

4VHe 

c  v,. 


n/v. 


Cv2-(0/p)2]  eizpe  dp 


p[p2e2i6-du2/c2+K+(p)][p2e2i6-o>2/c2+K  ( p)  3 


(A3) 


2  2 

First  note  that  neglecting  the  ca  /c  terms  in  the  denominator  in 

2 

Eq.  (46)  is  equivalent  to  neglecting  (v^/c)  compared  to  unity.  Then 
the  magnetic  field  from  the  branch  cut  is  given  by 


B  .ii 

Br  io>  dz 


1  r  °*P  -2ip 
-  r  ojM-  J  “re 
4  o  2 
c 


I 


[v2-(n/p)2]  eizpe  dp 

0/vf  CP2+K+e-2i6l[p2-fre-2i6] 


(44) 


To  evaluate  the  fields  at  the  surface,  substitute  u  *  ft/v^p 
and  z  -  0  in  Eq.  (43)  and  Eq.  (44)  to  yield 


2  2  i 

■J  “  mn  T  1 

EBr<°>  ■•JV^V  3U 


lo 


2  3 

(1-u  )u  du 


[!  +  (£)  kV^x  ♦  (£)  5V“] 


(45) 


3  m  -216  * 

v°>  -  • !  -rf * 

c  n 


J. 


(1-U2)u2  du 


[1  ♦  (V)  ♦  (TS*)  -  ^ 

v  2 

But  |(^)  K*|  is  0(A),  and  therefore,  to  lowest  order  in  A 


(46) 


EBr<°>  '  •ll“oJ1nie'3‘6 


(47) 


BBr<°)  “ 


1  T  A  -2i6 
T  u  J  A  e 
^  o 


(48) 


The  magnitude  of  the  ratio  (R  )  of  the  contribution  to  the  electric 

o 

field  from  the  branch  cut  to  that  from  the  pole  at  the  surface  of  the 
metal,  to  lowest  order  in  A,  is 


Using  Eq.  (30)  to  express  the  ratio  as  a  function  of  <ju,  the  maximum 

value  of  the  above  expression  can  be  shown  to  occur  at  m  ■  x>! 

2  2 

(provided  v  «  which  is  well  satisfied  for  most  metals).  The 

3 

maximum  value  of  R  is  «  .03  (<d  v£/uc)  .  For  silver,  (R  )  ■ 

opr  o  max 

13 

R  (3.6  x  10  )  «  .06.  For  poorer  conductors,  (R  )  will  be 

o  r  *  x  o  max 

considerably  less.  The  relative  contribution  from  the  branch  cut  also 

falls  off  rapidly  with  frequencies  far  from  Again,  for  light 

at  10^  cm  *  shining  on  silver,  Rq(2  X  10 mt  .005. 

A  complete  evaluation  of  E  (z)  requires  integration  of  Eq.  (43) , 

Br 

which  can  only  be  performed  numerically.  However,  it  is  possible  to 
find  an  asymptotic  representation  for  Egr(z)  valid  for  (Oz/v^)  »  1 
by  successively  integrating  by  parts  (v^/0  <  v^/v;  for  silver  this 
representation  is  always  valid  for  z  ~  10  ^  cm). 

,,  V,  ,,  h  iZa5^Vf  -zu/v, 

EBr«-T-oJ*f r-‘  < 


The  ratio  of  the  electric  fields  in  the  metal,  R  .  to  lowest  order  in 


A,  is  given  by 


(z) 


„„  ,  .  ■z',/v£  In(kc)z 
24R  (co)e  e 


Ep(z) 
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provided  Oz/v^  »  1.  It  is  interesting  to  observe  that  for  optical 

frequencies,  Im(k  )  ~  <u  /c  *  .5  X  10^  cm  and  with  exception  of  the 
c  p 

best  conductors,  (tu  v f/vc)  <  1.  Hence,  for  z  sufficiently  large,  R  (<u) 
p  t  z 

can  be  greater  than  unity.  However,  this  occurs  when  both  |E^|  and 

|E,  I  are  negligibly  small.  For  <ju/v  «  1,  R  (tu)  «  1,  and  the  electric 
1  b  1  z 

field  falls  off  exponentially  with  z. 


THE  REFLECTION  COEFFICIENT 

The  analysis  has  shown  that  when  an  electromagnetic  wave  is 
incident  upon  a  metal,  energy  is  transmitted  in  the  metal  as  a  com¬ 
plicated  electromagnetic  disturbance.  The  contribution  from  the  branch 
cut  is  dissipated  within  a  distance  of  the  order  v^/fi.  The  contribu¬ 
tion  from  the  pole  has  a  longer  range,  and  its  behavior  is  essentially 
determined  by  the  classical  model,  with  a  small  correction  (above 
liquid  helium  temperatures)  due  to  the  nonzero  Fermi  velocity.  The 
fraction  of  energy  reflected  from  the  surface  can  be  calculated  by 
applying  the  usual  boundary  conditions  at  the  metal  surface: 


EJ_  +  Er  -  E(0)  -  Ep(0)  +  Eg(0) 


(52a) 


Ei  -  Er  -  cB(0)  -  cBp(0)  +  cBg(O)  (52b) 

where  E^  and  Ef  are  the  incident  and  reflected  electric  fields, 
respectively.  The  constant  J  can  now  be  computed  by  adding  Eqs.  (52a) 
and  (52b): 

4E 

j  -  -  (i + f  [(* + r>lco + 2e"216li  <53> 
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As  A  -»  0,  Eqs.  (35)  and  (53)  yield 


lyo>l 


(54) 


the  classical  result  for  the  amplitude  of 
metal  surface.  I 


The  reflection  coefficient,  r 


the  transmitted  field  at  the 

2 

,  to  lowest  order  in  A,  is 


given  by 


r_ 

r 


1  + 


2Ltnt  [■ 


2  2 
(l-nr)  +  n^ 


2  2 
(1-n  y  +  nT 
'  r'  i 


-] 


(55) 


where  r  • 
c 


1-n 


1+n 


,  0  «  -[sin(25)  +  %sin  a>],  and  nf  and  n^  are  the 


real  and  imaginary  parts  of  n  ,  respectively. 

c 

2  2  2  2  2  *5 

For  a>  »  as  »  u  ,  the  correction  term  is  0(v,/c  )  ~  10  . 

P  i 

2  2  2 
05  »  1)  »  05 

P 

r/r  ~l  +  4-^ 

C  05p 


For 


(56) 


which  again  gives  a  negligible  correction.  Reuter  and  Sondheimer,^ 
using  a  similar  analysis,  have  discussed  the  reflectivity  of  silver 
at  liquid  helium  temperatures  (corresponding  to  much  larger  A's  in 
this  analysis)  and  have  concluded  that  the  correction  to  the  classical 
result  would  still  be  very  difficult  to  measure. 
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IV.  SECOND  HARMONIC  GENERATION 

In  this  section,  the  formalism  developed  in  Section  III  is 

2 

utilised  with  the  new  assumption  that  terms  of  order  are  no  longer 
negligible.  The  electromagnetic  fields  are  expanded  in  a  Fourier 
series. 


£-  l 


E  e 
— n 


-irvnt 


B  -  S'  B  e 

L 


-ina>t 


(57) 


n»-co 


and  solutions  are  sought  for  which  the  series  converge  rapidly.  These 
representations  are  valid  when  E^,  oscillates  with  frequency  o>  and 
has  a  sufficiently  small  amplitude.  In  this  case  each  of  the  nth- 
order  Fourier  coefficients  contains  a  term  which  is  the  nth  power  of 
the  incident  wave,  and  the  series  converges  quickly  for  small  amplitude 
waves. 

Substituting  Eqs.  (2)  and  (57)  into  Eq.  (1),  the  second-order 
component  of  the  transport  equation  becomes 

df„  _  3f  .  df. 


9  o  ft  tk  V  1 

(v-2io)f2 


(58) 


Let  of  * 


of  *  cd  +  %i\>  and  define  g^  by 


m  ao 

!i  *  L  dvx  L  %  v* 


(59) 


Then  multiplying  Eq.  (58)  by  v^  and  Integrating  yields 


-2  W 


dg.  •  •  df 

*1  +  v«  sT  ‘  »  E2J  J.„  J  Jvy  *1  5^  -  “l 


(60) 
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where 


•  *  3f.  3f,  Sf1 

D.  *  —  f  f  \E.  - —  +  B.  fv  - —  -  v  - —  )  ]  v,  dv  dv  (61) 

1  m  L  1  1  V  x  3vv  z  5v  /  J  1  y  x 


and  summation  over  j  is  implied.  From  the  relation  b fQ/3Vj  *  ( v j/v) 
(dfQ/dv),  it  follows  that  the  third  term  in  Eq.  (60)  vanishes  unless 
i  ■  j.  Further,  since  there  is  no  second  harmonic  component  in  the 
incident  wave,  the  boundary  conditions  (Eq.  (55))  require  that  the 
second  harmonic  electric  field  in  the  metal  be  longitudinal,  or  j  ■  z. 
For  the  z  component  of  Eq.  (60),  the  third  term  integrates  to 


3 

2 


(62) 


Similarly,  f^  (given  by  Eq.  (5))  is  an  even  function  of  v^  and 

an  odd  function  of  v  ,  from  which  it  can  be  shown  that  D  =  D  =  0, 

x  x  y 

and  D  reduces  to 


z 


eBj  ® 

- v 

m  z 


w  w 

u. 


v  r —  dv  dv 
x  Sv  x  y 

no  z 


(63) 


The  equation  for  g  seems  intractable  as  it  now  stands  since  the  exact 
z 

forms  of  the  fields  and  distribution  function  of  the  fundamental  are 
complicated  by  boundary  effects.  Fortunately  it  is  possible  to  make 
reasonable  approximations  which  will  render  Eq.  (60)  tractable.  In 
the  last  section  it  was  shown  that  the  branch  cut  contribution  to 
the  fundamental  field  has  a  small  amplitude  and  short  range,  so  only 
the  contribution  from  the  pole  will  be  included.  A  second  assumption, 
which  is  implicit  in  the  classical  theory  and  its  conclusions,  is  that 
the  electric  field  may  be  regarded  as  spatial?.y  constant  for  the 


V 


purpose  of  calculating  the  current  at  a  point.  This  assumption  is 
valid  when  the  field  does  not  change  appreciably  over  a  distance  of 
the  order  of  the  mean  free  path  of  the  conduction  electrons,  that  is 
to  say,  when  the  mean  free  path  is  small  compared  to  the  effective 
field  penetration  depth.  This  condition  is  equivalent  to  setting 
A  =  0  and  is  satisfactory  for  all  frequencies  of  interest  except  at 
very  low  temperatures.  The  solution  of  Eq.  (4)  can  therefore  be 


written 


which  yields 


ieE.(z)  3f 
1  _ o 

TOo'  dV 


AsiE 

mm'  1 


<z>  J 


3f  , 

v  - —  d  v 
x  9vx  - 


for  the  current  density.  When  Eq.  (64)  is  substituted  into  Maxwell's 

equations,  the  classical  dispersion  relation,  Eq.  (30),  is  recovered. 

Then,  using  Eqs.  (33)  and  (64)  and  carrying  out  the  indicated 

integration  results  in  the  following  equation  for  g  : 

z 


where 


88s  21*'*  .  ... 

ai - v  lz  •  v2 

z 


iV 

2  3 
mv^ 


[F(z)-E2(z)“| 


F(Z)  "  E^(z)B^(z) 


The  solution  of  Eq.  (66)  that  satisfies  the  specular  reflection 
condition  at  z  ■  0  is 


tCVo  l,,  ' 


for  z  <  0  and 


+ 

g2 


v 

z 


2io£/v 
z  z 
e 


e 


-2ioj|/vz 


21 <u|/v 

z 


(69b) 


for  z  >  0. 

Following  the  same  procedure  used  to  derive  Eq.  (13)  yields  for  the 
second  harmonic  current  density 


h 


05 

^(z+D  *(l)d|  -  J  ^(z-D  *U)d|] 


(70) 


where 


2  v 
3io>:  .  f 


GQV , 


Jo 


2U&"  |xj/vz 

e  v  dv 

z  z 


(71) 


With  this  definition  of  K2(x)  ,  and  introducing  Maxwell's  equation  for 
the  longitudinal  harmonic 


E2  '  21o)C. 


hmQ 


(72) 


results  in 


2mv"p.  ®  ® 

E2(z)  +  ('3neXJ0K2(Z'5)  i(5)d5  *  j  KOd|]  -  0  (73) 


or,  using  Eq.  (68), 


00  CD 

E2(z)  -  Q*  ^(z-l)  E2(5)d|  -  j^z+J)  E2(|)dtJ 


W(z)  (74) 


F 
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where 


W(z)  -  J”  ^(z-O  F(-|)d|  -  1^(2-!)  F(5)d|  (75) 


Since  W(z)  ■  -W(-z),  E2(z)  is  antisynmetric,  and  Eq.  (74)  becomes 


E2(z)  -  f  KjCz-l)  E2(|)d5  -  W(z) 

*7  -<T5 


(76) 


Fourier  transforming  Eq.  (76)  yields 


-m  Wfk)  Cf+c«  -  f+'k>l 

E  (k)  -  — “ —  -  lL(k)  - - - 

l-K^k)  1  l-K^k) 


(77) 


where 


^(k)  -  J  ^(x) 


2  2 
lo  u"  v,  v  dv 
-ikx .  p  r  f  z  z 

e  dx  -  -  ^  J 


and 


v^oc  k2  JO  v2  -  (2co'/k)2 


F.(k)  -  f  F(z)  e”ikzdz 
J0 


(78) 


(79) 


As  in  the  linear  case,  the  branch  cuts  in  the  k  plane  are  those  regions 

where  K^k)  is  not  defined,  i.e.,  where  kvf/2a/  is  real  and 
kvf  ~  1 

1 2^|  >  1.  K^k)  is  analytic  in  the  rest  of  the  plane  and  is  given 
by 


~  r  ,1  +  kvf/203# 

K2(k)  ■  ■  v2<uk2  [l  •  kv£  lnU  -  kvf/2oj')i 


where  In  denotes  that  branch  of  the  logarithm  which  vanishes  for 
k  -  0. 


(80) 


tf 

! 

i 


From  Eqs.  (68)  and  (79) 


where 


F.(-k)  -  F  (k)  -  a  -5 - j 

+  +  k2-(2kc)2 


«  -  -  -*7  E? 

1  c(l+nc)2 


Referring  to  Eq.  (79),  it  can  readily  be  seen  that  ^(k)  has  no  isolated 

singularities  in  the  cut  plane,  and  that  the  singularities  of  E^Ck) 

are  therefore  the  poles  of  F  (-k)  -  F  (k)  and  the  zeros  of  1-K_(k). 

■f  £ 

The  expression  for  F  (-k)  -  F  (k)  indicates  that  there  is  a  pole  of 

+  4* 

E„(k)  which  yields  a  signal  with  wave  number  k  =  2k  which  is  twice  that 

JL  C 

of  the  fundamental  frequency,  and  since  it  is  a  second  harmonic  term, 
it  therefore  has  the  same  phase  velocity  as  that  of  the  fundamental. 

The  roots  of  1  -  ^(k)  are  discussed  in  Appendix  B.  In  the 
limit  CD  »  u,  it  is  shown  that  there  always  exists  one  and  only  one 
zero  of  this  expression. 

Finally,  in  Appendix  C,  the  relative  amplitudes  of  the  residues 
are  calculated,  and  the  amplitude  of  the  branch  integral  is  estimated 
at  z  -  0.  It  is  shown  that  the  amplitudes  of  the  waves  from  the 
residues  are  nearly  equal  at  z  ■  0  but  are  approximately  180°  out 
of  phase.  When  cu  «  cup  the  amplitudes  at  z  ■  0  are  «  a  and  -or  from 
the  residues  at  1  -  K^k)  »  0  and  k  -  2kc,  respectively.  While 
both  values  of  k  correspond  to  evanescent  waves  for  v  «  ou  «  n>p,  the 
latter  has  a  much  longer  range  (~  c/wp)  than  the  former  (~  vf/oup)  and 
is  the  only  disturbance  of  interest  away  from  the  surface.  In  this 
frequency  range 
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M 


4ieE, 


mco3 


1.67  X  10'13  Ej  (mks  units) 


(83) 


for  a)  »  u.  The  propagating  wave  with  k  =  2k  is  smaller  in  magnitude 

^e^i  -10 

than  the  fundamental  in  the  ratio  ss - r  («  10  E.  for  o>  in  the 

mcoy  l 

infrared)  and  has  half  the  range.  The  heating  effect  of  the  second 

harmonic  is  therefore  down  by  approximately  the  square  of  this  ratio. 

When  the  inequality  o)  «  0)  is  not  satisfied  (e.g. ,  03/03  ^  1/3)  the 

P  P 

expressions  for  a  and  for  the  residues  are  more  complicated  and  the 
expressions  derived  in  Appendices  B  and  C  must  be  used. 
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V.  OBLIQUE  INCIDENCE 

For  the  case  where  the  propagation  vector  of  the  incident  wave 
makes  an  angle  0  with  respect  to  the  normal  to  the  metal  surface 
(i.e.,  the  z  direction),  the  harmonic  of  the  distribution  function 
and  the  components  of  the  fields  are  proportional  to 

-iqajz  sin  8/c  -iajt 
e  e 


Assuming  specular  reflection  and  following  the  procedure  outlined 
in  Section  111,  the  first-order  distribution  function  is  found  to  be 


df  Vvz  z  -ia*/v 

*+  _  e _ o  e  f  Tf  ^  *  E 


ft 4 - (v  ff 

1  mv  9v  v  l  x  l_L 


»  Ui£/v 

xU>dl  +  ^  •  Ex(8)d|] 


+  V 
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z  -ico$/v 


iaj|/v 


+  JQ  e  Z  Ey(|)d5] 


a, 


z  -i <D|/v, 


1<d|/v 


+  v  |  I  e 

z  U0 


(|)d|  -  e  Z  Et(5)d|1}  (84a) 


for  v  >0  and 
z 


-  e  o  * 

f,  ■  -  —  - - 

1  mv  Sv  v_ 


j  ■ 

l 


-i4D5/v, 


Ek(l)d? 


k-1 


(84b) 


2  2  2  % 

for  <  0,  where  v  *  (v^  +  vy  +  vj  .  Using  Eq.  (10)  and  introducing 
Maxwell's  equations  yield  a  set  of  three  integral  equations  for  the 


components  of  the  electric  field  which  can  all  be  solved  by  Fourier 
transforms.  The  results  yield  propagating  components  of  the  fields, 
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which  are  given  by  Che  usual  classical  or  Fresnel  equations  with  small 
corrections  of  order  L,  plus  transient  components  of  short  range  and 
small  amplitude. 

In  discussing  the  second  harmonic,  the  spatial  derivatives  of 
f^  may  again  be  ignored  in  the  "source  tenra"  for  f2  which  yield  three 
additicnal  integral  equations  for  the  components  of  the  second 
harmonic  fields.  The  boundary  conditions  for  oblique  incidence  no 
longer  require  a  reflected  harmonic  to  vanish,  and  the  amplitude 
may,  in  principle,  be  calculated  by  the  method  described.  However, 
even  using  the  Fresnel  relations  for  the  fundamental  fields  yields 
extremely  tedious  first-order  differential  equations  for  the  second 
harmonic  current  generators.  Again,  in  order  to  match  boundary 
conditions,  it  is  necessary  to  finu  the  residues  at  z  *  0  from  the 
poles  of  the  transformed  kernels  as  well  as  those  at  twice  the 
fundamental  propagation  vector  components.  These  have  been  shown  to 
be  approximately  equal  in  magnitude.  (The  branch  cut  contributions 
will  again  be  small.) 

Using  Eq.  (83),  the  ratio  of  the  energy  flux  reflected  with 
frequency  2®  from  the  surface  to  the  incident  flux  can  be  estimated 
from  the  present  analysis  to  be 


h  _  /4eE 

P.  \mcrn 
1  p 


.  v  CD 

:)  F<9-  & 


(85) 


where  q  is  the  incidence  angle  and  F  is  of  order  unity  with  F(0,  o>  /a/)  * 

P 

F(n/2,  a)p/ oq 0  “  0.  For  an  incident  laser  beam  of  ~  2.5  x  lQ^v/meter, 
the  fraction  of  the  incident  power  reflected  in  the  second  har -.onic 
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is  ~  10  .  For  order-of-magnitude  estimates  of  P  /P.,  a  plot  of 

2 

(4eE./mcuu  )  versus  incident  electric  field  E.  is  shown  in  Fig.  3. 
i  p  i 

For  electric  fields  ^  10^  v/meter,  the  inetal  will  probably  be  damaged. 

(3) 

Brown,  Parks  and  Sleeper,  '  using  a  ruby  laser  with  1  MW  peak 

power  and  <  50  nsec  duration  shining  on  silver  have  obtained  reason- 

4 

able  agreement  with  F  «  cos  0.  They  point  out  that  these  observations 

are  consistent  with  a  second  harmonic  polarization  proportional  to 

E^(V-E^).  This  expression  has  arbitrarily  been  extracted  from  the 

known  nonlinear  polarization  for  a  free  electron  gas  where  boundary 

effects  are  not  included,  which  contains  an  additional  term  proportional 

to  E^  x  B  .  and  which  yields  a  much  more  complicated  dependence  on  0. 

(12) 

Jha  *  has  attempted  to  include  the  boundary  effects  for 
oblique  incidence  by  using  an  iteration  procedure,  perturbing  about 
the  classical  Fresnel  solutions.  However,  for  the  fundamental  wave, 
his  correction  term,  to  lowest  order  in  v^/c,  is  of  the  same  order  of 
magnitude  as  the  Fresnel  solution.  This  result  shows  that,  first  of 
all,  the  iteration  procedure  is  invalid  near  the  surface.  Secondly, 
the  results  of  the  analysis  presented  here  indicate  that  the  total 
correction  to  the  Fresnel  solutions  are  of  order  A,  and  one  should 
obtain  only  the  classical  result  as  v^/c  -*  0. 

For  the  second  harmonic,  Jha  attributes  the  E^(V»E^)  term  to  a 
surface  phenomenon,  the  E^  X  Bj  expression  to  a  volume  contribution, 
and  obtains  a  complicated  expression  for  F(0,  ui^/io')  in  the  absence 
of  collisions.  These  correspond  to  the  root  of  l-i^Ck)  =  0  and  the 
residue  at  k  =  2kc ,  respectively.  While  such  a  partition  enables  one 
to  estimate  the  correction  in  the  metal  fairly  well,  the  procedure 
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again  breaks  down  near  the  surface.  Here  Jha  has  neglected  terms  in 
his  iteration  which  are  of  the  same  order  of  magnitude  as  those  he 
has  retained.  Basically,  as  has  been  shown,  there  are  three  terms 
to  consider  for  the  second  harmonic.  While  the  branch  cut  contribu¬ 
tion  is  indeed  negligible  at  the  surface  for  w  «  to  ,  it  is  necessary 

P 

to  include  the  residues  from  the  transformed  kernels  in  order  to 
obtain  better  than  the  order-of-magnitude  agreement  predicted  by 
Eq.  (85). 

The  results  presented  are  sufficient  to  obtain  rough  estimates 
of  nonlinear  heating  and  harmonic  generation.  While  the  amplitudes 
are  small,  it  may  still  prove  necessary  to  account  for  the  second 
harmonic  generation  in  applications  where  instabilities  might  be 
excited  and  grow  to  such  amplitudes  as  to  affect  the  operation  of 
a  device.  The  utability  of  a  system  when  a  wave  of  frequency  2iu 
is  generated  should  therefore  be  considered  and,  when  necessary,  the 
growth  rate  should  be  estimated. 
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Appendix  A 

ROOTS  OF  THE  LINEAR  DISPERSION  RELATION 

With  the  restriction  that  f  be  a  symmetric  function  of  v^, 
a  valid  representation  for  K^,(k)  with  the  same  region  of  analyticity 
as  Eq.  (23b)  is  given  by 


(w '  2 

n  c  -® 
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r  r  %  r 


dv 


®  2  1  -  (kv  /u)')2 


The  zeros  of  D(k)  are  then  roots  of 


(A-l) 


i  -  -  rrr  — ^ ^ 

V»«  JJJ  !  .  (kv  >■ 


dv  dv 

y  z 


(A-2) 


,  /\2 


Let  (kvf/<u')  =  a  +  bi,  where  a  and  b  are  real,  and  define  p  =  v/u). 

2  2  2 

Then  (kc/ou)  =  (c/vf)  (a  +  bi)(l  -  p  +  2ip)  .  The  general  technique 

for  locating  the  zeros  of  D(k)  involves  calculating  the  change  in  the 

argument  of  (D'(k)/D(k))  around  a  suitably  chosen  contour  in  the 

k  plane.  However,  for  all  metals  except  at  liquid  helium  tempera- 

13 

tures,  p  is  small  for  frequencies  greater  than  ~  10  ops.  With  this 
restriction,  the  real  and  imaginary  parts  of  Eq.  (A-2),  for  p  =  0,  are 


a(c2/v2)  =  — ^  JJJ  f0(vx’V^  1" - 2^2^ - 2~{\ 

f  n  ^  JJJ  o  x  y  +  (b^2)2J 


dv  dv  d*i  (A- 3) 

x  y 


and 


-b(c2/v^)  =  -Ej  JJJ  !  (.,,.  #)  [• 

r»  rat  ' 


JL 


n  ou 
o 


(l-an)“  +  (bu) 


l]  dvxdvy^  (A"4) 


respectively.  The  integrand  in  Eq.  (A-4)  is  always  positive,  and  this 
equation  can  then  only  be  satisfied  if  l<  =  0.  Equation  (A-3)  becomes 
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i  (  2/  2. 

1  -  a(c  /vf)  -  —  ^ 
no“ 

which  can  readily  be  solved  by  sketching  the  left-  and  right-hand 

sides  of  the  equation  versus  a.  The  left  side  is  a  straight  line  with 

y-intercept  unity  and  a  negative  slope  of  (c/vf)  .  The  right-hand 

o 

side  is  zero  at  a  «  •«  and  has  y-intercept  (lu^/io)  .  Inspection  of  the 
denominator  of  Eq.  (A-5)  reveals  that  the  slope  is  always  positive  for 
a  <  1.  (The  region  a  >  1  corresponds  to  the  branch  cut.)  Hence  there 

exists  one  and  only  one  intersection  for  p  »  0. 

2  2 
For  (aip/o))  <  1,  the  root  a^  is  positive  with  0  <  a^  <  (v^/c)  «  1. 

When  co  <  o>p,  the  root  is  negative  and  Ja^  \  can  easily  be  seen  to  be  less 

13 

than  1  provided  (OpVf/coc  <  1,  or  for  frequencies  greater  than  ~  10  cps, 
which  was  our  original  assumption.  (This  is  a  sufficient,  but  not  a 
necessary  condition.) 

! 

For  p  nonzero  but  small,  the  root  will  have  a  small  imaginary 
i  part,  but  its  magnitude  will  remain  <  1. 


^ -  dv  dv  dp. 

au  x  y 
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Appendix  B 

THE  ZEROS  OF  K2(k)  =  1 
The  roots  of  1  -  ^(k)  =  0  are  the  zeros  of 


by  term  to  yield 


teas'  ^  -i 

1  2. 
p  x  dx 

(3- 1, 

2  J 

CD 

JO  kv,  2  . 

P 

1  ■  fe?)  * 

1  the  integrand  may  be  expanded  and 

integrated  term 

kv  2 

( — 2: 
\2afJ 

(B-2) 

u  CD  1 

P 

O  I 

The  root  given  by  Eq.  (B-2)  is  not  valid  when  Icd  cd"/cd  I  &  ttt  or  »  1. 

p  i  iu 

The  latter  case  is  of  little  interest.  A  solution  in  the  former  case 
may  be  found  by  noting  that  a  valid  representation  of  Eq.  (B-l)  is 


%s  to* 

2 

CD 


3p2[l*p  tan  1(l/p)l 


(B-3) 


2  2  - 1 

where  p  *  -  (2a//kvp  ,  and  that  branch  of  tan  (1/p)  is  chosen 

which  equals  n/2  at  p  ■  0.  Then  for  |cd  a/7a)p|  «  1,  Eq.  (B-3)  may 

be  expanded  about  p  »  0  to  yield 


(B-4) 


2 

Let  (kvf/2o/)  *  a  +  bi,  where  a  and  b  are  real,  and  consider  the 

collisionless  limit  y  «  0.  Then  of  “  co  and  equating  the  imaginary 
parts  yields  b  3  0.  Sketching  the  left  and  right  sides  of 
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£D_ 


2 . 

x  dx 


I* 

-ax 


(B-5) 


reveals  that  one  (and  only  one)  intersection  always  exists  whenever 
the  right-hand  side  of  Eq.  (B-5)  is  defined  (i.e.,  -®  <  a  <  1), 
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Appendix  C 

THE  RESIDUES  AND  BRANCH  CUT  CONTRIBUTION  FOR  THE  SECOND  HARMONIC 


The  residue  of 

2Qfk  K^k) 

[k2-(2kc)2][l-K2(k)] 


(C-l) 


•t  l-KjCk)  ■  0  is  given  by 
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First  consider  the  range  where  Eq.  (B-2)  is  valid.  Then 
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When  |ui  ao'/tu^l  «  1,  we  use  Eq.  (B-4)  which  gives 
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and  the  residue  is 
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The  residue  of  Eq.  (C- 1)  at  k  =  2k^  is 


aK^k.) 

1-K2(2kc) 

When  Eq,  (B-2)  is  valid  the  residue  is 
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where  s  »  ~  (-7—1  (1  -  — *^7)  can  be  neglected. 
5  \<a  c/  \  <0  cdV 

When  Eq.  (B-4)  is  valid  the  residue  is 
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In  a  manner  identical  to  that  used  to  derive  Eq.  (41)  ,  the 

contribution  to  the  branch  cut,  E„  ,  is 
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where 
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.  Then 
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With  the  substitution  u  =  2fJ,/v^p,  Eq.  (C-10)  becomes 
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have  assumed  »  1.  Then  Egr(0)  CWJfl'/cOp.  Hence  the  ratio  of  the 
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contributions  from  the  branch  cut  to  the  residue  is  ~u)ft'/(n  . 
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model  and  for  normal  incidence,  the  second 
harmonic  waves  in  the  metal  are  calculated 
and  the  magnitude  of  the  relative  heating 
is  computed.  The  method  for  solving  the 
general  problem  for  an  arbitrary  incidence 
angle  is  outlined  and  the  ratio  of  tie 
average  energy  flux  in  the  second  harmonic 
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to  the  incident  flux  is  estimated. 


